Model Independent Tests of Cosmic Growth vs Expansion 
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We use Gaussian Processes to map the expansion history of the universe in a model independent 
manner from the Union2.1 supernovae data and then apply these reconstructed results to solve 
for the growth history. By comparing this to BOSS and WiggleZ large scale structure data we 
examine whether growth is determined wholly by expansion, with the measured gravitational growth 
index testing gravity without assuming a model for dark energy. A further model independent 
analysis looks for redshift dependent deviations of growth from the general relativity solution without 
assuming the growth index form. Both approaches give results consistent with general relativity. 



I. INTRODUCTION 

Cosmological surveys have advanced to provide preci- 
sion determinations of the distancc-rcdshift relation and, 
to a lesser extent, the growth-redshift relation for cosmic 
history at redshifts z < 1. These improvements allow 
not only determination of cosmological parameters but 
fundamental tests of the cosmological framework in a 
more model independent manner. Rather than assum- 
ing a model with cold dark matter plus a cosmologi- 
cal constant (ACDM) or plus dark energy parametrized 
by a constant equation of state ratio w or time varying 
w(z) = wq + w a z/(l + z), one might like to investigate 
the expansion history H(z) and growth factor D(z) or 
the growth rate f(z) = d\n D /d\n(l + z) directly, with 
minimal assumptions. 

Reconstruction of the expansion history, in terms of 
the inverse Hubble parameter H~ 1 {z) or deceleration pa- 
rameter q(z) = — din H / dln(l + z) — 1, can be carried 
out purely kincmatically, without assuming a particular 
theory of gravity or field equations (i.e. Fricdmann equa- 
tions). Gaussian Processes [1] prove to be an effective 
statistical technique for carrying out such a reconstruc- 
tion from distance data, as done in [2]. 

Growth of matter density perturbations into large scale 
structure, however, depends explicitly on the dynamics, 
i.e. the gravitational force law. Within general relativ- 
ity (and pressureless matter being the only significantly 
clustering component), expansion and growth are locked 
together, either one determining the other. Given that 
recently growth data have advanced to cover a reason- 
able redshift range, z m — 0.8, at ~ 10% precision, it 
is interesting to test whether this interrelation actually 
holds. We can enlarge the Gaussian Process technique to 
do this in a model independent manner (rather than as- 
suming a dark energy parametrization), although some- 
what less generally than the previous expansion history 
reconstruction in that we must separate out the matter 
density. 

In Sec. II we briefly review the cosmographic recon- 
struction of expansion history from distance data and 



describe the extraction of growth history from large scale 
structure data, in particular using redshift space distor- 
tion measurements. Section III carries out likelihood 
analysis of current data and derives confidence contours 
on the matter density and gravitational growth index. 
We outline future applications and conclude in Sec. IV. 



II. GAUSSIAN PROCESS METHOD 
A. From Distances to Expansion 

Distance measurements play an essential role in our 
understanding of the history and contents of the uni- 
verse. They have a linear relation to (through integra- 
tion of) the inverse Hubble parameter or expansion rate 
for a spatially flat Roberston- Walker universe as we as- 
sume. This linearity is an important property for a Gaus- 
sian process (GP), since the derivative (or integral) of a 
GP is another GP, making error propagation particularly 
straightforward. Thus, for a luminosity distance 

d l (z) = (l + z) V (z) = (l + z) f dz'H-^z'), (1) 

Jo 

if we model di as a GP then the conformal distance r\ or 
inverse Hubble parameter H~ x is one as well. 

Gaussian processes provide a robust statistical method 
for using stochastic data measured at certain points (red- 
shifts) and reconstructing the full function (distance- 
redshift relation or inverse Hubble parameter) describing 
the underlying relation, complete with covariances and 
without assuming a specific model for the relation. Sec 
[1] for detailed explanation of their general application, 
and [2-5] for specific application to dark energy and cos- 
mology (also see [6, 7], though they fix several aspects of 
the GP and distance model, and [8] for a genetic algo- 
rithm approach). Here we follow most closely [2]. 

Given data y at a set of points Z we reconstruct the un- 
derlying function f , or its derivatives, at any set of points 
Z\ . The probability distribution functions arc Gaussians 
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described by a mean function m(Z) and covariancc ma- trix k(Zi,Zj): 
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The inferred mean and covariance of the functions are 
given by 



X^(Z,Z)y (4) 
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For the GP covariance function we use a common form, 
the squared exponential, 

fc(z, z') = o-j exp f- ^ V (6) 
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where cr / defines the overall amplitude of the correlation 
and I measures the coherence length of the correlation. 
The parameters a 2 and / are hyperparameters in the fit. 

If f(z) is the reconstructed distance then /' is the re- 
constructed inverse Hubble parameter. Within general 
relativity, H^ 1 also determines the linear growth his- 
tory of large scale structure. While the growth factor or 
growth rate will not be linear functions of H , and so 
are not GPs themselves, the error propagation is still di- 
rect. The basic approach is that supernova distance data 
allow (model independent) GP reconstruction of H~ 1 (z), 
with its covariances between redshifts, as in [2], and then 
this can be propagated to predictions of growth. These 
can then be compared to growth data from galaxy red- 
shift surveys. 



B. From Expansion to Growth 

The linear growth factor is difficult to measure directly, 
free from astrophysical effects such as galaxy bias. Weak 
gravitational lensing data, which docs not involve galaxy 
bias, is not currently sufficiently accurate to be useful 
for the desired reconstruction. Therefore we use galaxy 
rcdshift survey measurements of the growth rate through 
redshift space distortions, whose anisotropic angular de- 
pendence allows separation from galaxy bias. 

Redshift space distortions arise as follows. The matter 
density perturbations forming large scale structure in- 
duce gravitational potential inhomogeneities, and these 
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in turn give rise to motions of the matter, or peculiar ve- 
locities. These velocities add to the galaxy redshift due to 
cosmic expansion, causing an anisotropic observed den- 
sity field, in redshift space. Since the peculiar velocities 
are proportional to the growth rate / = dhiD/dhia, 
where the scale factor a = 1/(1 + z), then these rcdshift 
space distortions can be used as a cosmological probe 
[9]. In the linear perturbation limit, [10] showed the ob- 
served (redshift space) galaxy power spectrum is related 
to isotropic real space density power spectrum by 



P s (k,n) = (b + f f i 2 ) 2 P r (k) 



where k is the wavemode of the density perturbation, \x 
is the cosine of its angle with respect to the line of sight, 
and b is the galaxy bias. 

Since the power spectrum is proportional to the square 
of the mass fluctuation amplitude, P r (k,a) oc erf (a) oc 
D 2 (a), then the redshift space distortion observable is 
fas oc dD/d\na (see, e.g., [11]). Normalized to the 
present mass fluctuation amplitude os.cb an excellent ap- 
proximation to the cosmological and gravitational depen- 
dence of this quantity is 



0(a) 



/0"8 
08,0 



n m (^e^ dlM 'l n »(°') 1 - 1 ] , (8) 



where 7 is a constant called the gravitational growth in- 
dex [12]. For general relativity and ACDM, 7 = 0.55. 
The gravitational growth index form has been shown to 
be accurate at the 0.1% level for a wide variety of dark 
energy and gravity models [12, 13], so long as the gravita- 
tional strength remains scale independent and no strong 
clustering of dark energy occurs. 

This form immediately allows us to carry out a test 
of gravity without choosing a dark energy model or 
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parametrizing the Hubble expansion, since 



(9) 



The normalization relative to today, i.e. the cts.o m 4>i 
ensures that there is no dependence on H(z) for red- 
shifts higher than the highest growth measurements (and 
hence distance data, since these extend further). We can 
thus use the model independent GP reconstruction of 
H^ 1 {a)j Hq 1 from the SN distance data, propagate it to 
predictions of the growth relation, and by comparing to 
the growth data fit for the matter density today fi mj o 
and 7, the latter testing gravity. 

In a second approach, one can actually enhance the 
model independence by writing 



cj)(a) = G R.(a) + 8cj)(a) , 



(10) 



where 4>gr fixes 7 = 0.55 as from general relativity, and 
use Gaussian processes to reconstruct the function 5(f) 
without assuming a functional form given by the growth 
index 7. (Mathematically, one uses the GR relation as 
the mean function in the GP and sees if the hyperparam- 
eter tr/ giving the amplitude of deviations is consistent 
with zero.) This allows for exploration of a wider variety 
of extended gravity theories. We use both approaches in 
the next section. 



III. RESULTS FROM DISTANCE AND 
GROWTH 

We apply GP to the Union2.1 compilation of super- 
nova distance data [14] and incorporate the /<rs(z) data 
from the Baryon Oscillation Spectroscopic Survey (BOSS 
[15]), SDSS DR7 [16], WiggleZ [17], 2dF [18], and 6dF 
[19] galaxy surveys. Note one must be careful to use the 
growth data values derived without assuming a specific 
expansion model. Through a scan over the likelihood 
surface, marginalizing over hyperparamctcrs, we can de- 
rive confidence contours for £l m fi and 7, or study the 
hyperparamctcrs themselves. 

Figure 1 shows the joint 2D contours on f2 mj o~7 for two 
different values of er^o- The results are consistent with 
GR value of 7 = 0.55, and the variation of o^o slides 
the contours in f2 m! o with little effect on the probability 
distribution of 7. 

The GP method, without assuming any dark energy 
model, indicates that 7 can take values in a considerable 
range, though general relativity i.e. 7 = 0.55, is right 
near the peak of the likelihood. In particular, even when 
fixing il m ,a = 0.28, say, values of 7 as low as those char- 
acteristic of scalar-tensor theories such as f(R) gravity 
[20] or as high as that of DGP gravity [21], 0.42 and 0.68 
respectively, are allowed at 95% CL. 

Current data therefore does not have the leverage to 
look for more subtle redshift dependent deviations from 
GR that might not be captured by 7, at least not with- 
out assuming a specific model. We quantify the reach 




FIG. 1. 68% and 95% joint confidence limits on 7 and il m ,o 
are shown derived without assuming a dark energy model, us- 
ing current supernovae distance and galaxy clustering growth 
data. The left contour of each pair has o^o = 0.801, the 
WMAP7 concordance value [22], and the right has crg.o = 0.78 
to show the effect of a small shift. 



of current data in two ways. Figure 2 shows the GP 
reconstructions of the growth rate as a function of red- 
shift, with data from BOSS, SDSS DR7, WiggleZ, 2dF, 
and 6dF ovcrplottcd. The light green band is composed 
of samples of reconstructions with A\ 2 < 3 relative to 
the best fit, when fixing 7 = 0.55, while the dark red 
band allows 7 to float. In both these cases we have fixed 
cr 8 ,o = 0.801. 
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FIG. 2. Reconstruction of the growth rate fas(z) is shown 
for the case when fixing 7 = 0.55 (light green curves) or al- 
lowing it to float (dark red curves). Current growth data is 
overplotted and we fixed <7s,o = 0.801. 

The second method involves taking the even more 
model independent approach of fitting for an arbi- 
trary time dependent correction to the general relativity 
growth rate, S(j>(a) — 4>{a) — 4>gr{cl)- That is, we take 
GR to provide the mean function for the GP and let the 
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FIG. 3. Limits on the amplitude of deviation a/ from the 
GR growth relation (with model independent dark energy) 
are shown in 2D joint confidence contours with the present 
matter density. 



data constrain the amplitude of the deviations given by 
the hyperparameter <jf. Figure 3 shows the 2D bound in 
the <7/-f2 mj o plane. The data prefers no deviation from 
general relativity, i.e. Of = is within 68% CL. 



technique for such model independent analyses. The GP 
reconstruction of the expansion history was juxtaposed 
with growth rate data from redshift space distortion mea- 
surements in galaxy surveys to obtain probability distri- 
butions involving the gravitational growth index 7. The 
general relativity value was found to be a good fit, al- 
though due to uncertainty in the matter density TOi o 
and to a lesser extent the mass fluctuation amplitude 
(78,0 a wide range of values is tenable within current con- 
straints. 

We further extend the model independence by look- 
ing for any deviation in the growth rate, without using 
the growth index formalism. Building on the GP recon- 
struction we test the growth data for deviations from the 
prediction of general relativity as a function of redshift. 
The results are again consistent with standard gravity, 
tested without assuming any particular model of dark 
energy. Stringent exploration of the laws of gravity, how- 
ever, requires more accurate future growth and distance 
data. Ongoing and future surveys will greatly enhance 
our ability to carry out model independent investigation 
of the cosmological framework. 
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IV. CONCLUSIONS 

We have demonstrated a method to solve for the ex- 
pansion and growth histories of the universe simultane- 
ously, without assuming any model or parametrization 
for dark energy. This is a key test of the cosmological 
framework since within Einstein gravity one determines 
the other. Using the results we derive from Type la su- 
pernovae and large scale structure data, we test for and 
quantify deviations from general relativity in two ways. 

Gaussian Processes (GP) provide a useful statistical 
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